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ABSTRACT 

Binary systems subject to generic perturbations evolve on quasiperiodic or- 
bits. We derive the most generic class of perturbations, which allow to evaluate 
secular effects via generalized complex true and eccentric anomaly parameters, 



^ ■ by use of the residue theorem. Such perturbations include both the generic 



Brumberg force and various linear contributions to the conservative dynamics of 



r-- 

\Q • compact binaries, up to the second post-Newtonian order. As an example, we 



compute the self-spin contribution to the luminosity due to gravitational radia- 



m 
o . 

tion of a compact binary consisting of galactic black holes, the most important 
type of source for LISA. This translates to simply calculate the residue of the 
instantaneous energy loss in the origin of the complex parameter plane, which 
q . illustrates the power of the presented method. 

Subject headings: compact binaries, perturbed Keplerian motion, parametriza- 



tion 



1. Introduction 



The perturbed two-body problem is one of the central areas of interest of Celestial 
Mechanics. A wide class of perturbations of the dynamics of a binary system are described by 
the generic Brumberg force [Brumberg (1991) ]. Another type of perturbation is provided by 
general relativistic modifications of gravitational dynamics. Such measurable effects arise in 
planetary motion (like the general relativistic contribution to the precession of the perihelion 
of the planet Mercury), however more important are the modifications in the strong- field 
regime, typically appearing in compact binaries composed of neutron stars and / or black 
holes. A well-known example is the Hulse- Taylor double pulsar [Hulse & Taylor (1975)], 
the orbital period of which is continuously decreasing as gravitational radiation escapes from 
the system [Taylor, Fowler, McCulloch (1979)] (for a recent review see [Weisberg & Taylor 
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(2004)]). Efforts to directly detect gravitational waves with Earth-based interferometric 
detectors (LIGO[Abramovici et al. (1992)], VIRGO [Bradaschia et al. (1990)], GEO [Hough 
(1992)] and TAMA [Kuroda et al. (1997)]) are currently under way. Gravitational waves 
from distant coalescing galactic black hole binaries will be searched for by the forthcoming 
Laser Interferometer Space Antenna (LISA) [Bender et al. (1996)]. With improving detector 
sensitivity, upper limits for the signal emitted by different configurations of gravitational 
wave sources were already set ([Abbott et al. (2004)], [Messaritaki (2005)], [Abbott et 
al (2005)]). Capturing gravitational waves in the noisy background requires a 3.5 post- 
Newtonian (PN) accuracy both in the orbital phase [Blanchet (2002)] and in the energy loss 
[Poisson (1995)] occurring in a compact binary. 

The backreaction of the escaping gravitational radiation on the orbit decreases the 
orbital period of the compact binary. Although higher order corrections to this effect are too 
small to be directly detected in the early stages of the inspiral, they are important in the 
latter stages, when the post-Newtonian parameter increases up to 0.1 (at 10 gravitational 
radii separation) and higher order contributions become quite significant. Therefore the 
contributions to the energy loss < dE/dt > and angular momentum loss < dL/dt > have 
been worked out up to 2PN order accuracy in the wave generation (or equivalently, to 4.5 
PN order following the Keplerian picture), in the case of eccentric orbits, and up to even 
higher orders for circular orbits. (Here <> denotes the time-average over one period, e.g. a 
secular contribution.) The corrections arise at each PN order from both general relativistic 
modifications of the dynamics and/or as a manifestation of various physical characteristics 
of the system, like the spins, mass quadrupole momenta and magnetic dipole momenta of 
the constituents. The final aim of gravitational wave astronomy is to incorporate all of these 
physical characteristics into accurate templates to be used in data processing. Recovering 
the physical characteristics of astronomical objects from a simultaneous analysis of data 
collected in both the electromagnetic spectrum and gravitational wave bandwidth would be 
a major breakthrough in modern astronomy. 

Driven by this motivation, the interest in the perturbed two-body problem has been 
renewed. In this context, in [Gergely, Perjes, Vasuth (2000)] a generic perturbed two-body 
problem with radial equation 

i=0 ^ 

r 2 N = 2- + 2 — 2 

fx r ji 2 r 2 

was considered. The point mass \i (representing the reduced mass of the binary) is orbiting 
a fixed gravitational centre m (the total mass of the system) at distance r (the separation of 
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the components). The energy E and the magnitude of the orbital angular momentum L are 
constants of the perturbed motion 1 . The expression f 2 N (the subscript N denoting Newtonian) 
has the familiary Newtonian functional form, however E and L characterize the perturbed 
motion. The perturbing terms contain the small coefficients tpi, assumed constants. As will 
be shown in detail in Section 3, both the generic perturbing Brumberg force, the first post- 
Newtonian order (1PN) general relativistic corrections [Damour & Deruelle (1985)] and the 
spin-orbit (SO) interaction for compact binaries [Barker & O'Connell (1979)] fit into this 
scheme. The generic equation (1) defines a radial period T and turning points r max . 

min 

With these assumptions the integrals of the type 

I(u,n) = £^dt (3) 

(with constant uj and n arbitrary integer) could be easily evaluated by use of the residue 
theorem. Important steps in the proof of [Gergely, Perjes, Vasuth (2000)] were the intro- 
duction of suitable generalized true and eccentric anomaly parametrizations of the radial 
motion, r (%) and r (£) as 

2 1 + cos x 1 — cos x , , n 

- = - + - , (4) 

T '"min '"max 

2r = (1 + COsO^min + (1 - COsO^max , (5) 

together with their complex counterpart z = exp (ix) and w = exp(i£), respectively. The 
advantage of this method over other methods for computing the secular integrals (3) is 
its striking simplicity. As proven in [Gergely, Perjes, Vasuth (2000)], in the majority of 
physically relevant cases the only pole is in the origin. Extremely rarely, a second pole (also 
given) occurs. The method was successfully applied in [Gergely, Perjes, Vasuth (1998)] 
for the evaluation of the averaged energy and angular momentum losses due to the SO 
interaction for compact binaries on elliptic orbit, and the results were in agreement with 
those of [Rieth & Schafer (1997)], obtained by a different method. 

The more accurate study of the inspiral phase of compact binaries however reveals 
other important contributions to be taken into account. These are the perturbations due to 
the (a) spin-spin (SS) interaction [Barker & O'Connell (1979)]; (b) the interaction of the 
magnetic dipole moments for neutron stars with huge magnetic fields (DD) [Ioka & Taniguchi 
(2000)] and (c) the contribution arising from the motion of one of the binary components 
(treated as a monopole) in the quadrupolar field of the other (QM) [Poisson (1998)]. The 



1 Separating E and L from the rest of the coefficients of r° and r 2 in the radial equation is advantageous, 
as it will enable us to express all results in terms of constants with physical interpretation. 
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corresponding secular contributions to the energy and angular momentum losses could be 
computed [Gergely (2000)], [Vasuth, Keresztes, Miliary, Gergely (2003)], and [Gergely & 
Keresztes (2003)] by the method presented in [Gergely, Perjes, Vasuth (2000)], in spite of 
the radial equation being different from Eq. (1), which rendered these computations outside 
the domain of validity of Theorems 1 and 2 of [Gergely, Perjes, Vasuth (2000)]. Therefore 
the question arose, whether the results presented in [Gergely, Perjes, Vasuth (2000)] would 
hold more generically. 

In all above mentioned cases the radial equation was a generalization of Eq. (1) with 
the coefficients tpi being periodic functions of the true anomaly x- Thus in this paper we 
raise the question, whether the results of [Gergely, Perjes, Vasuth (2000)] would hold for 
generic periodic functions ipi (x) • The answer we find is that some restrictions apply, and we 
derive the class of admissible periodic functions ipi. Section 2 contains the announcement 
and proof of the generalized Theorems 1 and 2, as compared to [Gergely, Perjes, Vasuth 
(2000)]. 

Starting from the explicit expressions for the respective interactions, we show in Section 

3, that the SS, QM and DD perturbations fit into the class of admissible perturbing functions 

Vi- 
lli Section 4, as a simple application of the presented method, we give the self-interaction 

spin contributions to the luminosity of a compact binary consisting of black holes, in the 

case when S 2 < 5i. 



2. Admissible perturbations and the integration method 

The turning points of the radial motion (1), with ifii = ipi (x) are the loci f = 0, given 

by: 

i-2 

(6) 

Here A represents the magnitude of the Laplace-Runge-Lenz vector characterizing a Kep- 
lerian motion with E and L and we have introduced the notations 

V?r = ¥>i(0), = • (7) 

Unless in the case presented in [Gergely, Perjes, Vasuth (2000)], here ip~ ^ (pf. The 
expression of the turning points allow to write up the parametrizations (4) and (5) in detail, 
whenever needed. 
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The integrals (3) could be formally rewritten in terms of the integration variable x as 

f 2 * u ( 1 dt \ , 



2.1. The case n > 

For n > the factor r ra in the denominator of the integrand in Eq. ( 8) has the binomial 
expansion 



2 



r(x) 



\ ' mm / \ ' max / 

By applying the chain rule dt/dx = r^dr/dx, with dr/dx derived from the true anomaly 
parametrization (4), and inserting r" 1 as the series expansion of the radial equation (1), we 
find for the second factor of the integrand (8) 

1 dt n f L 2 \ 

^dx' L\ + 4^A 2 sm 2 X ) ' (10) 

r = ±U + Ai mx -*em , (id 

i=0 ^ ' 

±<pt(Gmn-Aoy] . (12) 

The second term of the expression (10) contains sin 2 x i n the denominator, and apparently 
becomes singular at x = kit (k any integer). The key point of the proof of Theorem 1 of 
[Gergely, Perjes, Vasiith (2000)] was that a factor of sin 2 x could be also separated from T, 
so that the integrand of / (u, n) became regular everywhere. The proof has also showed that 
there is no other x-dependence left in the denominator. 

We would like to see whether this convenient property of T holds for generic periodic 
functions ifi (x). For this we decompose ifii ( X ) as follows 

oo 

^ (x) = + 9a cos X) sinJ X , (13) 

3=0 

the coefficients and being constants. Then the condition 



k / \ n— k 



<pf = fio ± gm 



(14) 
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follows. Eq. (13) is equivalent with a decomposition into a generic Fourier series, however 
better suited for our purposes. Its advantage lies in the property of the terms with j > 2, of 
manifestly containing the factor sin 2 x- Therefore singular integrands in (8) could arise only 
from the terms containing f i0 , g i0 , fn and gn. We will discuss these one by one. 

As fio is the constant part of (pi (%), the contributions containing f i0 are the constant 
perturbing functions already dealt with in [Gergely, Perjes, Vasiith (2000)]. Thus the terms 
containing f i0 are regular. 

A similar algebra as for the terms with f i0 applies to the terms containing g i0 . By 
employing the binomial expansion (9) and the expression of the turning points (6), those 
contributions to T which contain the factor 2 (fi/L 2 ) 1 g iQ sum up to 



£({) (Gm^A 



l-(-l)'... 



sin x + 2/+i(x) 



(15) 



The expressions E k (x) are defined by Eqs. (29)-(31) of [Gergely, Perjes, Vasiith (2000)], 
and they are proportional to sin 2 x- This proves that the terms with g i0 do not give singular 
contributions to the integrand of / (u>, n). 

What remains to check are the j = 1 terms of T, given by 

^ 2(/ii + &icosx) sinx ^ 
i=o r 

(We have used that, according to Eq. (14), the expressions A± receive no contribution from 
fn and gn.) As the coefficients fn, gn are already of first order, the Keplerian true anomaly 
parametrization r = (yu/L 2 ) (Gm/x + Aocosx) can be employed in Eq. (16). A tedious but 
straightforward algebra, based on (a) the binomial expansion of l/r l ; (b) the replacement 
cos 2 x = 1 — sin 2 X-, whenever possible; (c) a second binomial expansion of (1 — sin 2 x) 9 , 
whenever applicable and (d) grouping together the terms without a sin 2 \ factor, into terms 
proportional to sin \ and sin x cos x, leads to the following two conditions on the coefficients 
fn and g a : 

J2 JT % (Gmfi ± A o y (f n ± 9il ) = , (17) 

i=0 

which have to be satisfied in order dt/dx to be regular for every x- 

The last question to be addressed is what are the conditions to be imposed on u, in 
order that I (u,n) stays regular. The answer is simple: oo can be any periodic function of Xi 
which is regular. Such functions will not affect the proof presented above. 
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With these conditions satisfied, all factors in the integrand of I (cu, n) become simple 
polynomials in sinx = [z 1 — 1) /2iz and cosx = [z 2 + 1) /2z, provided n > 0. This renders 
the only pole in the origin, and we have proven: 

Theorem 1: For all perturbed Keplerian motions characterized by the radial equation 
(1), with periodic perturbing functions ipi obeying the conditions (13) and (17), and for 
arbitrary periodic functions u the integral I (cu, n > 0) is given by the residue in the origin 
of the complex true-anomaly parameter plane. 



2.2. The case n < 

In order to cover the complementary case of n < 0, we employ the eccentric anomaly 
parametrization (5). The integrals I (u,n) then are evaluated as 

'<".»>= r^Cz^u. us) 



l r n+1 \rd£ 

and for any n! = — n — 1 > we apply the binomial expansion 

n' 

(20 B ' = £ ft) ^nnC'ax^l + COsO*(l " COsO"'"* , (19) 
fc=0 

which is a polynomial in cos£. From the eccentric anomaly parametrization (5) and the 
radial equation (1) (with <p (x (£)) i n place of <p (%)) we find 



t = £ fnt-rt-«»{-^' ) • (2D 



i=0 



Qi = (-^j^t^^-Ao)- 2 

^-(Gm/z+Ao)-! . (22) 

Here (p~ and </?^~ are identical with those defined in terms of the true anomaly (this is, because 
at r m ax both angles x and £ take the values and 7r, respectively). Considerations completely 

min 

analogous to those in the proof of Theorem 1, together with the relation between the two 
parametrizations 

Gmfi cos ^ - A 

cos X = ^ ~ A 7, (23) 

Gm/j, — A cost, 

(-2£L 2 //i) 1/2 sin£ 
Gmfi -A cos£ v ; 
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after lengthy computations lead to the following result. Provided the Eqs. (17) hold, a factor 
of sin 2 £ separates out from T and cancels out the sin 2 £ factor of the denominator in Eq. 



As main difference with respect to the true anomaly parametrization, the terms ifi (£) /r l ~ 
in the expression (21 ) could lead to additional factors (Gm/j, — A cos£) in the denominator 
of the integrand. Therefore, the second pole 



could appear in the process of passing to a complex eccentric anomaly parameter. We can 
now enounce our 

Theorem 2: For all perturbed Keplerian motions characterized by the radial equation 
(1), with periodic perturbing functions ifi obeying the conditions (13) and (17), and for 
arbitrary periodic functions u>, the integral I (u,n < 0) is given by the sum of the residues 
in the origin of the complex eccentric- anomaly parameter plane and (for any with i+j > 2 
or gij with i + j > 1 in (fi) at Wi . 

Our Theorems contain as special cases the Theorems presented in [Gergely, Perjes, 
Vasuth (2000)]. There, the conditions (17) are automatically satisfied, as (fi are constants. 
It is remarkable, that the constancy of y?j can be relaxed in a quite generic way, in which the 
infinite number of coefficients in the expansion of tpi (%) should obey only the two contraints 



In this Section we present two important classes of dynamical systems, for which our 
Theorems apply. 

3.1. Perturbations characterized by the generic Brumberg force 

The generic perturbing Brumberg force [Brumberg (1991)], [Soffel (1988)], includes a 
wide range of perturbations employed in Celestial Mechanics. It can be derived from the 
Lagrangian 



(20). 




1/2 



(25) 



(17). 



3. Admissible dynamical systems 



Cn + &PB 



(26) 
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£w = ^ + G^_ (27) 
If o , ^\ 4 Gmji 2 



A\ Gmu 9 / A\ G 2 m 2 a 



where a, /3, 7 and A are perturbation parameters. The radial equation derived in [Gergely, 
Perjes, Vasiith (2000)] from this Lagrangian is of the form (1), with the small coefficients 

c 2 ^ = _3(2a - 2(3 + \)E 2 , 
cVf = -4Gm/i(3a - 2/? + 2\)E , 
cVf = -2 / u 2 G 2 m 2 (3a-2/3 + 2A + 7) 

+ — (2a — 2/3 + X)EL 2 , 
/i 

cVf = 2Gm(a + 2A)L 2 . (29) 

As all ipf are constants, these perturbations fit into the domain of validiy of the methods 
discussed in [Gergely, Perjes, Vasiith (2000)]. 



3.2. Orbital evolution of compact binaries 

Another important example is provided by compact binary systems, consisting of black 
holes / neutron stars. Their dynamics is conservative up to the 2PN order, the first dissipative 
effects due to gravitational radiation occuring at 2.5 PN order. The conservative dynamics 
including all leading order contributions, due to either various physical characteristics, or to 
the general relativistic modifications of the motion is described by the Lagrangian 



■>CB 



C N + C PN + C S o + £ss + £qm + £dd , (30) 



with the various contributions derived first in [Damour & Deruelle (1985)] (PN), [Keresztes, 
Mikoczi, Gergely (2005)], (SO), [Kidder, Will, Wiseman (1993)] (SS), [Poisson (1998)] 
(QM) and [Ioka & Taniguchi (2000)] (DD) read 2 : 



r 1 n a \ A, Gm ^ 



(3+rj) v 2 +r]f 2 



2 The magnitude and direction of the spins are denoted as Si and Si. The angle subtended by them 
is 7 = cos _1 (Si • S 2 ). The total spin is S = Si + S2 and er = (m 2 /mi)Si + (mi/m 2 )S2- The mag- 
nitude and direction of the magnetic dipolc moments di are denoted as di and dY They subtend the 
angle A = cos _1 (di • d 2 ) with each other. In a coordinate systems /C with the axes (c,L x c, L), where 
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Cso = ^jv • [ r x (4S + 3(7)] , 

(Si • S 2 ) - ^ (r • Sx) (r • S 2 ) , 

£dd = ^[3(n-d 1 )(n-d 2 )-d 1 -d 2 ] . (31) 

Note that the 1PN accurate relativistic corrections to the Keplerian motion character- 
ized by Cpn is the particular case of the perturbing Brumberg Lagrangian Cpb, with the 
specifications a = rj/2, j3 = (1 + 3r))/2, 7 = 2 + 77, A = 2 — 77 (where 77 = \ijm). Therefore the 
PN perturbations fit into the domain of validiy of the methods discussed in [Gergely, Perjes, 
Vasuth (2000)]. 

The same holds for the SO perturbations, where the only non-vanishing tpf° in the 
radial equation is 

so = JfV ( 4L . S + 3 L . CT ) . (32) 
cr 

This is a constant, as the angles /% span by the spins with the orbital angular momentum L 
are constants (to leading order) 3 . 

The situation is entirely different for the SS, QM and DD contributions to the dynamics 
of a compact binary. In all these cases the energy E is a constant of motion, while (in contrast 
with the 1PN and SO contributions) the magnitude L of the orbital angular momentum is 
not. Indeed, with the notable exception of the PN contribution, it can be derived from the 
respective Lagrangians that the orbital angular momentum L = r x p is not conserved: L 7^ 



£ss — 



G 

c 2 r 3 



-QM 



Gfim" 
2r 5 



^2 Pi 



c is the unit vector in the J x L direction, the polar angles Ki and tpi of the spins are defined as 
Si =(sinKiCOSipi,smKiSmipi,cosKi) (see [Gergely, Perjes, Vasuth (1998)]). In the coordinate system A? 
with the axes (bi, Si x b;, Si), where bi are the unit vectors in the Si x L directions, respectively, the polar 
angles on and j3i of the the magnetic dipole moments di are di =(sin on cos j3i, sin on sin cos on) (see [Vasuth, 
Keresztes, Mihaly, Gergely (2003)]). The quadrupolar parameters (see [Gergely & Keresztes (2003)]) are 
defined as p, = Qi/m^m 2 , where Qi is the quadrupole-moment scalar [Poisson (1998)] of the i th axially 
symmetric binary component with symmetry axis Si . 

3 Note however that neither £50 nor the radial equation with the perturbing function (32) characterize 
uniquely the SO-perturbation. This is, because a spin-supplementary condition (SSC) should be additionally 
imposed. We have followed here the SSC of [Pryce (1948)] and [Newton & Wigner (1949)]. There are other 
known SSC-s, including the covariant SSC employed in [Kidder (1995)], which gives the non- vanishing 
constant perturbations ip^°' cov = (2E/c 2 m) (L-er) and ipf°- cov = - {2G^/c 2 ) (2L • S + L-cr). These can 
again be expressed in terms of the angles Ki , which are constants in the leading order, required here. 
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0. While for SO perturbations this means merely a precessional change about the vector 
4S + 3-cr (so that L is conserved), the situation is more complicated in the SS, QM and 
DD cases, in which L is a function of the orbital position L = L(x). In all these cases 
however an angular average L over the radial motion can be computed, which is conserved 
to the accuracy of linear effects we consider here. The reason for this is that up to 2PN, 
the dynamics is conservative. A description in terms of the constants of motion E and L 
is then in order, as described in detail in ([Gergely (2000) ], [Gergely & Keresztes (2003)] 
and [Vasiith, Keresztes, Mihaly, Gergely (2003)]. The radial equation in all these cases can 
be expressed in the form (1), with L in place of L while the perturbing terms contain the 
following explicit x-dependences: 

^ss _ _J^_ g 1 g 2 s [ nKl smn 2 {2A cos (x+ 8)] 
2c 2 L 

+ (3Gm/2 + 2Acosx)cos(2x+ S)} , 
G 

(fin S = -^ r S , iS , 2 f3cos/ti cos/t 2 — cos 7 
& 

—3 sin K\ sin k 2 cos(2x+ 8)} , 

<P 2 M = - Gm _^ V V sin 2 Ki{2 A cos(x+ Si) 
2L 



QM 
f3 



+ (3Gjj,m + 2Acosx) cos(2\; + <5;)} , 
= —Gfi 2 m 3 '^2Pi 1 — 3sin 2 /tj cos 2 (x+ 



rf D = -^^[(SGm^ + AAcosx)B 2 (x)-AsmxB' 2 ( X )}, 
L 

y?f D = -pd x d 2 [A - 3B 2 (x)} , (33) 

where ipo is the angle between the periastron and the node line [Gergely, Perjes, Vasiith 
(1998)], ip = (ipi + ip 2 ) /2 represents an average, 8 = 2<^ — tp) and 8i =ipQ — ipi differences 
in the azimuthal angles. The quantty A is the magnitude of the Laplace-Runge-Lenz vector 
characterizing a Keplerian motion with E and L. Finally, A® and B 2 (x) denote 

A = 2 cos A + 3(pia 2 - p 2 a l ) sin (8i - 8 2 ) 

-3(pip 2 + o"iO"2) cos (5i - 8 2 ) , (34) 

B 2 (x) = {vi°2 - P\Pi) cos(2x + 8) - (pi<j 2 + p 2 ax) sin(2x + 8) , (35) 

where pi = sin ctj cos Oi = sin oii sin fa cos /tj + cos «j sin Kj. Note that for the SS, QM, and 
DD interactions there are no contributions fi for any i ^ 2, 3. 

It is straightforward to transform the expressions (33) into the form (13), and to identify 
the coefficients fn and g n . From among the various emerging contributions, those entering 



- 12 - 



in the constraint (17) are / 21 = Aj, f 31 — , g 2 \ = Gmfi^y and gr 3 i = — ^L 2 /yuJ 7, with 7 a 
complicated expression depending on the particulars of the SS, DD or QM contributions 



7 

IQM 

iss 

1DD 



ISS + IQM + 1DD , 

6G/i 3 m 3 



sin 2 Kj sin 25j , 



6Gfi 2 SiS 2 ,. 

= n — sin Ki sin /t 2 sin , 

c 2 L 



(0-1 o- 2 - pip 2 ) sin (5i + 5 2 ) 



+ (PlC2 + P2C1) cos (5i + 5 2 ) 



(36) 
(37) 

(38) 



(39) 



It is easy to verify that the condition (17) holds in all cases, irrespective of the particular 
form of 7. This is why the method developed originally for constant yvs could be employed 
for these cases as well 



4. The self-spin contribution to the luminosity of compact binaries 

Gravitational radiation produces dissipative effects in the orbital evolution. While the 
energy E, and the total orbital momentum J are conserved up to 2PN order accuracy 
(including the conservative perturbations considered in this paper and the 2PN relativistic 
perturbation), they start to evolve due to escaping gravitational radiation. This dissipative 
evolution can be conveniently computed with our method, presented in Section 2. The 
secular losses due to gravitational radiation are exactly integrals of the type (3), in which 
the condition of periodicity of uj is obeyed for all losses of SS, DD and QM type. As already 
argued in Section 3.b, the radial equation written in terms of L is also of the form (1) in 
all of these cases. Therefore our Theorems can be applied in computing secular radiative 
changes. Such secular effects were discussed in ([Gergely (2000)], [Vasuth, Keresztes, Miliary, 
Gergely (2003)] and [Gergely & Keresztes (2003)]) for the SS, DD and QM contributions, 
respectively. 

In this Section we would like to apply the method for computing the so-called self- 
spin contribution to the gravitational luminosity (energy loss) of compact binaries composed 
of black holes. For a black hole the quadrupole moment scalar is directly related to the 
spin: Qi = —Sf/rrii [Poisson (1998)]. Therefore the QM losses are effectively self-spin 
contributions and add to those discussed in [Gergely (2000)], [Mikoczi, Vasuth, Gergely 
(2005)]. Such contributions become important whenever the mass ratio of the compact 
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Table 1: The coefficients in the self-spin contribution to the luminosity. 








1 


2 


3 





-539 784 


-1229 200 


-657120 


-47040 


i 


797076 


1807400 


959 280 


67680 


2 


266 615 


818 986 


601812 


60 600 



binary is at least rj = 0.1, which is frequently the case for colliding galactic black holes, 
which are the most important source types to be searched for by LISA. In such cases, as 
S 2 /Si oc if, we can safely ignore S 2 . Therefore the only SS interaction is the self-spin 
contribution. 

The instantaneous energy loss can be found by taking the appropriate derivatives of 
the mass and current quadrupole tensors of the mass //, as described for example in [Kidder 
(1995)]. The self-spin contribution to the instantaneous energy loss computed in this way 
is a complicated expression of r, f and x- By employing the parametrization (4), we obtain 
dE/dt as function of x alone. The secular contribution to the energy loss arises by integration 
of dE/dt over one radial orbit. This integration is carried out by computing the residues 
enclosed in the circle ( = e %x . Due to Theorem 1, we are assured that there is only one pole, 
at ( = 0. The summed-up self-spin contribution to the luminosity is 

G 2 {-2E/j,) 3 / 2 S 2 m 2 



-QM+SS 



960c 7 L n 



mi 



Z + 



Z x + Z 2 cos 2 



where 



'13 



3=0 



EL 2 
G 2 m 2 \i 3 



(i/>o-i/> 



0..2 



sin 2 k 



(40) 



(41) 



The angles k and ip are shown in Figl. The coefficients are enlisted in table 1. Eqs. 
(40)-(41) contain the totality of second order contributions in the spins to the luminosity, 
provided S 2 = 0, and as such, represent the first correction to the Lense-Thirring approach, 
discussed in [Gergely, Perjes, Vasuth (1998)]. 



5. Concluding Remarks 



Perturbed Keplerian motions can be parametrized in many different ways. (For a re- 
view of these parametrizations see [Klioner & Kopeikin (1994)].) The method of com- 





Fig. 1. — The figure represents a black hole - black hole binary system. S is the spin of the 
supermassive black hole with mass mi, while the spin of rri2 <C vn\ is negligible. As there is 
either no SS or QM-contribution to the orbital angular momentum [Kidder (1995), Gergely 
& Keresztes (2003)], L = /xr x v is the Newtonian orbital angular momentum, k is the 
angle between S and L. The angle between the node line and projection of the spin in the 
plane of the orbit is denoted as if). 

puting secular effects presented in [Gergely, Perjes, Vasuth (2000)] relies on novel complex 
parametrizations of the radial part of the motion and its advantage over previous approaches 
[Ryan (1996)], [Rieth & Schafer (1997)] is in its overwhelming simplicity. Other methods 
can of course be successfully applied, like the use of the Laplace second integrals for the 
Legendre polynomials ([Whittaker & Watson (1996)], [Gopakumar & Iyer (1997)]). 

In the present work we have generalized both the class of integrands and the class 
of perturbing forces for which the simple method relying on the use of complex true and 
eccentric anomaly parametrizations applies. 

In a generic perturbed two-body problem, it is easy to see whether a particular secular 
effect to be computed is of the type (3), with uj periodic, and then to check whether the 
periodic perturbing functions (fi obey the conditions (17). Then, depending on the power n, 
either Theorems 1 or 2 can be applied. 



We have also shown explicitly, that for several physically interesting linear perturbations 
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in the conservative dynamics these Theorems apply. Whether the SO, SS, QM or DD 
perturbations are of 1PN order or higher, depends on the specific system considered. For 
comparable mass compact binaries for example the SO contribution is of 1.5PN order, while 
the SS, QM, DD contributions are all of order 2PN. 

By contrast, when one of the masses of the compact binary dominates over the other, 
the SS contribution becomes negligible, however the so-called self-spin contribution of the 
dominant spin becomes important. Driven by this remark, as an application to the presented 
method, we have computed the self-spin contribution to the gravitational luninosity. 
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